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ABSTRACT 

in this paper we consider the problem of the joint design 
of multiple transmitters and receivers in a MIMO- 
OFDM sceoario. Wc consider that both the transmitters 
and receivers have multiple antennas* coofiguring a 
Multi-input-Muld-Ou^ut (MIMO) channeL Based on 
diis structure, space diversity techniques can be used to 
separate the signals transmitted by difTerent users. 
Concretely, the proposed approach consists in 
beamfonning at both sides of the system, combined 
with the Orthogonal-Ffequcncy-Division-Multiplexing 
(OFDM) modulation. The goal is the minimization of 
the total transmitted power subject to constraints related 
to the Quality of Service (QoS) for each user in terms of 
the Bit EiTor Rate (BER) and the maximum transmitted 
power for some terminals. By making use of this 
^proacb, the Multiple Access Interference (MAI) is 
minimized and the physical resources are exploited in 
an optimum way. 'fhe design is based on Simulated 
Annealing (SA), an heiirislic''6&hni<pic*3BrrT6"1!nd bie 
optimum design in tlus ttolihC» 
can take into accoimt atrkind of coiuitram^ " ' 

I. INTRODUCTION 

In last years mobile communications have experimented 
an exponentia] growth in tenns of capaci^ requirements 
or number of users, bit-rate increase and Quality of 
Service (QoS) or Bit Error Rate (BHR) inqirovement 
Space diversity based structures represent a good 
collection of techniques to cope widi these challenges. 
Concretely, the Multi-Input-Multi-Output (MIMO) 
channels, based on the utilization of multiple antennas 
at both the transnutters and receivers, have been the 
target of study of many researchers. Moreover, many 
new stancferds such as the European Wireless Local 
Area NetwoHc (WLAN) HIPERLAN/2 II] have 
proposed to use the Oitbogonal-Frcqucncy-Division- 
Multiplexing (OFDM) as then* modulation formaL 
Many works have considered the combination of MIMO 
and OFDM from a singl&-user point of view [2] [3], but 
this structure is also able to implement the Space 



Division Multiple Access (SDMA), in which several 
users accede to the radio channel at the same time and at 
the same frcqucacy, providtitg a very high spectral 
efficiency able to accommodate important increases in 
system capacity. 

In this work we propose a multiple mssst system, in 
which several Mobile Terminals (MTs) with more than 
one antoma are coexisting in the same area. These users 
are assumed to use the same frequency and accede to 
the channel at the sanoe time by means of the OFDM 
modulation. The receiver equipments^ for example die 
Base Stations (BSs) in an upliiik transmission, are also 
assumed to have more than one antenna and must detect 
the transmitted symbols by all the users. The approach 
taken into account in this work for designing the 
transmitters and receivers is joint beamforming. For 
example, in an uplink commimication with one BS, the 
transmitters use a bcamformcr for each frequency or 
carrier, and the receiver has a bank of beamformers for 
each carrier in order to detect the transmitted symbols 
by all the users. Because of simplicity and 
computational load reasons, we consider that linear 
singtc-user detectors arc implemented at the receiver 
side. 

The stated problem consists in the design of the 
beamformers subject to different constraints related to 
the Quality of Service {(}oS) for eacb communication in 
terms of the Bit Error Rate (BER). Our goal is to 
minimize the total transmitted power so as to guarantee 
an optimum utilization of the physical resources. As in 
the uplink channel some of the MTs may have a limited 
battery capacity, we also permit to add constraints 
corresponding to the maximum transmit power for some 
of the MTs. This optimization problem is not easy to 
solve and closed expressions do not exist Besides, it is 
not convex, aiKl so classical approaches based on 
gradient algorithms may find a suboptimal design. In 
this work we propose the ^plication of Sinmlatcd 
Annealing (SA) to solve this optimization probtem! 
13 an heuristic technique able to find the global bjgjinmm 
subject to all kind of ccmstiaints. 

This paper is organized as follows. In Section II 
the system and signal models arc presented. Section III 
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aiqslies SA to the stated design problem, whercas in 
Sectioo rv ao altemativs Lagrangogradicnt based 
algorithm is proposed Finally, in Section V some 
simulation results and conclusions arc obtained. 

IK SYSTEM AND SIGNAL MODELS 

A. System and Signal Models 

We consider a wireless scenario with several users 
acceding to the radio channel at the same time and the 
same frequency. A typical exan^le corresponds to an 
uplink transmission in a cellular system such as that 
shown in Fig, 1. 




Fig. 1. Typical configuration in a multiuser MIMO- 
OFDM scenario with 3 users or comnrnnicationSw 

The most general model would correspond to an ad- 
hoc network, m which all the devices can behave as a 
transmitter and receiver simultaneously, and each 
communication or link is engaged to a different 
transmitter and receiver, incMing possible multihop 
routing. In this paper we concentrate on the case of an 
uplink transmission in which K MTs arc transmitting 
sinniltancously in time and frequency to the same BS, 
that is responsible for detecting the iaroimation symbols 
transmitted by all the users. As it was stated before, 
we assume that both the MTs and BS aie equipped with 
more than one antenna. The extension to other kinds of 
scenarios or networks is direct and sunplc [4]. 

The snapshot vector signal model y^/) at the nth 
carrier or frequency at the BS is [2]: 

y»(0=tHl"biV(0+ii„(/) (J) 

where the number of elements of y^/) is equal to the 
number of antennas of the BS. b«*'* is the transmit 
beamvector corresponding to the /th MT and /rth carrier^ 
and has as many components as number of antennas in 
the associated MT. s}\() is the transmitted symbol by 
the Ah user at the wth cairicr during the fth OFDM 
symbol period. It is assumed that the information data 

are normali2xxl; ^ll^-i'^wfj^ I , where E\ \ stands 

for the mathcnstical expectation operator. The matrix 
H,*'' represents the MIMO channel between the lih MT 
and the BS at the /rth carrier p]. Finally, djc snapshot 
vector B„Cr) models (he contribution of noise plus 
interferences from outside the system at the BS at the 
same frequency. The associated covariancc marrix is: 



V„=^{n„(On^(/)}, where {f stands for 
conjugate transpose. 

B. Single-User Receiver Optimization 

The BS receives the signals transmitted by all the MTs. 
Its goal is to detect or estimate the information symbols 
transmitted by all the users. Naturally, the optimum 
receiver should be based on a multiuser detector; 
however, as its computalJonal load is too high and 
possibly prohibitive, in this work wc propose the 
utilization of a bank of linear single-user reccivcis, 
where each of them considers the other users as 
interferences. The optimum receive beamvector uj^^ for 
the ^h user at the nth frequency is the one that 
maximizes the Signal-to-Noisc aiKi Interference Ratio 
(SNIRX and corresponds to the well-known Matched 
Filter (MF) [2]: 

ar*),«i*)R(*r>„mb.M (2) 
RL**=9.+iHlX^rHi'>" (3) 

where a,^*> is a constant that does not affect the SNIR 
and that can be calculated so as to have well adjusted 
dedsk)n threshold levels. 

The estimates of the transmitted symbols are 

cafculated as follows: ii**(0 = <fec|al*^''y„(o}. 

where dec{^) stands for decision. Taking into account 
this design, the SNIR at the output of the beamformer at 
the /rth carrier for the kxh user can be expressed as: 

SNIlC' = bi*'"Hi*'"Rl,*>"'Hl'>bl*' (4) 
In OFDM it is defined the effective BER^'^ for the 
Aih user as the BER averaged over the subcairiere (here 
we assume that N carries are active for each user): 

where Ai„ is a constant depending on the modulation 
format for each subcarricr (for BPSK, Kh =2) and wc 
have assumed that the interferences are approximately 
Gaussian distributed. 

IlL SIMULATED ANNEALING BASED 
TRANSMFfTER OPTfMIZATION 

In last section we derived die optimum linear single- 
user receivers. In an uplink scenario, the BS has a bank 
of /T single-user receivers, each of them responsible for 
detecting the symbols of one user, while considering the 
others as interferences. In thfe section our goal is to 
design the transmit bcamivectors bw^*^ subject tn QoS 
constraints. In this paper these QoS are directly 
expressed in terms of the target effective BER. 
Obviously, an important parameter in a wireless 
scenario is the transmitted power, specially in an iq>link 
transmission or an ad-hoc network, in wWdi the 
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transmitters are MTs with a limited battery cqxacity. 
Because of this, wc also take into account these 
considerations in fite design. 

The power used for transmitting the signal at the 

mh frequency by the Ath user is proponional to |l>I,**[f • 

The optimization problem consists in minimizin^s Ihe 
lotaJ transmitted power Pj (6) subject to the QoS 
constraims, which Indicate that the effective BER for 
each user must be lower than or equal to a target 
prefixed value y*** (7), guaranteeing a minimum QoS: 

Pr=±%n (6) 

BEF^'^^y''' Jt = l„..,A: (7) 
In addition to the QoS constraints (7) and as 
explained before, it could be convenicm to add 
constraints corresponding to the maximum transmit 
powers for some of the MTs P„J\ specially in the 
uplink channel and when the batteiy capacity is veiy 
limited. Let us define T as the set of MTs to which the 
maximum transmit power constraints are ^plicd. This 
constraints are expressed as shown in (8): 

-0 (8) 

Closed expressions do not exist for this 
constrained optimization problem, and so, it must be 
solved by means of iterative methods. 

This problem formulation generalizes that 
presented in [5]. There, the design was based on a 
multiuser Multicarrier CDMA (MC-CDMA) system 
widi only one antenna at both the transmitter and 
receiver, and the QoS constraints were applied directly 
in tcnns of the SNIR. The problem was solved by 
means of a gradient-based algorithm. In [6] the 
transnaitted power, which was a prefixed value, was the 
same for all the users and carriers; therefore no power 
allocation was carried out The goal was the 
optimization of the SNIR by means of an iterative 
algorithm {Alternate & Maximize - AM). [7] 
generalized the iterative technique presented in [6] to 
more general conditions and a5$unq>tions. The main 
problem of [5J, [6] and [7] is that a non-optimum 
solatioii, (bat is, a local minimum may be achieved due 
to the non-convexity behavior of this optimization 
problem [7]. Besides, neither of them considered the 
individual transmit powers constraints (S), The gradient 
technique presented in [5] was based on the Lagrange 
multipliers method and the quadratic penalty function 
[8]. This technique may have some convergence 
problems related to the speed of convergence and the 
convergence itself [5] [8]. Besides* the gradient requires 
that the consn:aints are difTerentiablc, limiting its 
application. In the simulation results section, we 
conqiare a gradient-based technique witii the proposed 
al^rithra in this pq>er. based on Simulated Annealing 
(SA). 



In this paper we propose the application of SA to 
solve the stated optimization problem taking into 
account all the constraints previously presented. SA is 
an iterative algorithin able to find the global o^HiTinm 
TOlu fabh,' even w^ieri" the* probleih ' "is ncuvxonvex;' 
4^i;ovtcfed a feasible solutioii exists, fliat is, a solution 
that satisfies all the constraints simultaneously. If ' a 
Teasible solution does not exist, then the algdfitfam will 
not converge to any acceptable design. SA has analogies 
with the annealing in physics as explained in [9J. In our 
problem, in each step there is a collection of 

beamvectors which is caUed the current 

solution. Given the current solution (equivalent to a 
concrete particles arrangement or state in physics), a 
new solution or collection of beamvectors is proposed. 
If it is *i>cttcr" than the previous one, then it is accepted 
as the current solution. On the contrary, if it is "Vorse" 
then the proposed solution is accepted with a certain 
probability. This mechanism, called *'hill-climbing", is 
important as it avoids finding a suboptimal solution or 
local mininuun. The paramtAer that controls this 
acc^tancc probability is the temperature as in the 
case of anoeating in physics. The higher the 
temperature, the higiher the acceptance probability. The 
temperature T is lowered step by step, so that 
asymptoticaJly only 'Tjetter" solutions are accepted and 
a minimum is achieved. The meaning of "better" and 
*%rarse" is related to the definition of a cost function J( ) 
that depends on the transmit beamvectors (9), and which 
corresponds to the energy of a state in the annealing in 
physics. 

Here wc summarize the basic ideas of the SA 
algorithm that we prtyposc to solve Ac stated 
optimization problem. In this case, we s^ply the 
algorithm to a continuous solution ^ace. Besides, the 
cost fimctionX ) depends on the temperature 7(9): 

• Cost Junction definition: 

where (x)* = max {0,x}. 

• Proposed solution generation: 
bj,*>=bi*'+wl*> «=0,....Ar-l.A = I K 

M,'„''~CN{0,oil) (10) 

• Probability of occ^tance the proposed 
solution {Metropolis criterion) [9]; 

" (f(m-m'})) 



(9) 



Pr = cxp 



(U) 



System canting: the temperamre is lowered 
with an exponential profile: 

T^fiT fi=^0.99 (12) 
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Initially the lemperatuix; T must be high enough so 
that most of the proposed sohitions are acxcpted (in our 
case, we increase the temperature until the acceptance 
ratio is 95%). In this work we run 100 iterations pa- 
each value of 7*. The cost fimction is equal to ihc total 
transmit power plus a quadratic penalty term that takes 
into account if the efifectivc BERs are greater than the 
lequired ones, and if the individual transmit powers arc 
greater than those specified. As 7 is lowered, the 
penalty term is increased, and so, we asymptotically 
avoid solutions that do not fulfill the constraints. We 
make relative comparisons of the effective BER and 
transmit powers with the required values by means of 
the log(-) function, as e)(|>erimcntally we have observed 
that it behaves better than absolute comparisons. 
Nevertheless, other kind of comparison functions could 
have been used. The proposed sohitions are generated 
by applying con^ilex circularly symmetric Gaussian 
noise to the components of the bcamvectors. The 
acceptance ratio is checked per each value of T, In case 
it is lower than O.l for S limes, then the variance of the 
Gaussian noise is lowered by means of nn exponential 

rule {a\ *— O.QSCT^). This helps the algorithm to 

find with more precision the minimum as the 
temperature is lowered, and therefore, it increases 
significantly the convergence rate of the technique. 

IV. LAGRANGE-GRA0IENT BASED 
TRANSMITTER OPTIMIZATION 

In this section an alternative solution based on the 
Lagrange multiplier mediod and the quadratic penalty 
fimction [5] [8] is presented. For simplicity, here we 
only take into account die constraints corresponding to 
the BER. It can be easily shown that the minimum 
transmit power can be achieved when these constraints 
are fulfilled with equality, and so, we apply this result 
directly. The gradient algorithm is based on the 
definition of the following Lagrangian expression L\ 



log 



(13) 



An iterative process is applied so as to find the 
optimum weight vectors b„^*. The penalty parameter X 
is calculated by means of a gradient ascent technique, 
whereas the beamvectors are calculated with a descent 
gradient approach. In both cases the parameter step-size 
\l must be used: 

(M) 

bL*'<-bL*'-;iV^i (15) 

For space reasons, the expressions corresponding 
to V are not shown here, although they can be 

caaly obtained by using the results deduced io [5]. 

A pn^lera of this technique is that no "a-priori" 
value of |t can be obtained for assuring the convergence 
in a reasonable number of iterations^ and thai a local 



J 



minimum may be achieved iastead of the global one. 
Besides, as it is explained in [5] and [8], the speed of 
this gradient-based technique decreases tmportandy as 
the design approaches a solution fidfilHng the 
constiBints, and therefore, in a reasonable time, it cannot 
be assured that it is possible to fmd a design fulfilling 
the constraints with an allowable error. 

V. SIMULATION RESULTS AND CONCLUSIONS 

In this section we simulate an uplink channel with 5 
MTs and 1 BS. The OFDM modulation consists of 16 
canieis and both the BS and MTs have S antennas. The 
QoS constraints in terms of the mean uncodod effective 
BER are: 10 \ 1 0"^ and 1 0'^ for each user, and a=l 00. 

For all die figures presented in this section we 
show die behavior of die techniques versus the number 
of ''flops" in MATLAB, so that we can compare tiie 
algorithms directly on a computational load basis. 

In the first scenario we assume that the path loss is 
very similar for all the users. In Fig. 2 we show the 
power for the 3 users and the BER when no constraints 
on the individual transmit powers are applied. It is 
concluded that the SA algorithm is able to find a design 
fulfilling the QoS requirements. The power allocated to 
the first user is 8.45 W, and the total power is 20. 1 W. 




Fig. 2. Pcrformanoe of the Simulated Annealing 
algorithm in scenario I . 

As it can be seen in the figure above, during die 
first iterations the transmit powers allocated to the MTs 
are very high, and as a direct consequence, the effective 
BERs are extremely low. This solution is not efficient, 
as we can lower the transmit powers by increasing the 
BERs until the QoS constraints are fulfilled with 
equality. This behavior is shown in the figure. 

If we apply a power constraint to the first user 
equal to 8 W, then the results obtaiiked by SA are those 
shown in Fig, 3, We conclude that in this case the 
algorithm allocates 7.6 W to the first user, whereas the 
others increase their corresponding power consumption. 
The global transmit power has increased up to 20.8 W. 
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Fig. 3. Pcrfonnancc of the Simulated Annealing 
algofiehm in scenario I . Power constraint on MT I . 

In Fig. 4 and 5 we show the results for an scenario 
in which the third user has a path loss with respect to the 
first two users equal to 12 dB. Fig. 4 corresponds to the 
application of SA, whereas Fig. 5 corresponds to the 
Lograngian gradient-based algorithm with a \l parameter 
equal to 0.001. We conclude that with the same 
computational load, the SA can fuliin the constraints, 
whereas the gradient based technique decreases 
importantly the speed of convergence, as the solution is 
nearer from them. 




Fig. 4. Performance of the Simulated Annealing 
algoridim in scenario 2. 
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Fig, 5. Perfonnance of the gradient based algorithm in 
scenario 2. 

REFERENCES 

[1] ETSI, HTSi TS 101 47$ vlJ.I: Broadband Radio 
Access Networks (BRAN): HJPERLAN T)fpe 2 
(PHT) layer, ApriXlim. 

[2] A. Pascual lacrtc, A. I. Pdrcz-Ncird, M. A. 
Lagunas, "Jomt Beamfbnning Strategies in 
OFDM-MIMO Systems'*, in Proceedings IEEE 
ICASSP*02, Orlando (Florida), May 2002. 

[3] D. P. Palomar. Jotai M. Cioffi, M. A. Lagunas, A, 
Pascual lserte» "Convex Optimization Theory 
Applied to Joint Beamforming Design in 
Multtcarrier MIMO Channels", suhmitted to IEEE 
GLOBECOM*02. 

[4] A- Pascual Iserte, A. I. Pirez-Nena, M. A. 
Lagunas, "Joint Transceiver Optimization in 
Wireless Multiuser MIMO-OFDM Channels 
Based on Simulated Annealing". Invited paper, in 
Proceedings EUSIPCO'02. Tolouse. Sept. 2002. 

[5] Tat M. Lok, Tan F. Wong, "Transmitter and 
Receiver Optimization in Multicarricr CDMA 
Systems", IEEE Trans, on Communic. vol. 48, pp. 
1197-1207, July 2000. 

[6] K. K. Wong. R, S. K. Cheng, K. B. Letaief. R. D. 
Muich, "Adaptive Antennas at the Mobile and 
Base Stations in an OFDM/TDMA System", IEEE 
Trans, on Communic, vol. 49, pp. 195-206, Jan. 
2001. 

[7] J-H- Change L. Tassiulas, F. Rashid-Farrokhi* 
"Joint Transmitter Receiver Diversity for EfRcient 
Space Division Multiaccess", IEEE Trans, on 
Wireless Communic, vol. 1, no. ), pp. 16-27, Jan. 
2002. 

[8] Dimitri P. Bertsekas, Constrained Optimization 
and Laff-ange Multiplier Methods^ Computer 
Science and Applied Matfaematicsi, Academic 
Press, 1982. 

[9] P. J. M. van Laarhoven, E. R L. Aaits, Simulated 
Ameaiirtg: Theory and Applications^ Kluwer 
Academic Publish^ 1987. 



PAGE 34«4-RC«>AT9«3«00ai2:04:15P«. pastern D.yagMTUne,.8VR:08PTC«FXRF^-DH«:^^^^ 



